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On  the  horizontal  variability  of  the  upper  ocean 

Daniel  L.  Rudnick 

Scripps  Institution  of  Oceanography,  La  Jolla,  California 


Abstract.  The  last  decade  has  seen  a tremendous  increase  in  the  number  and 
quality  of  observations  of  horizontal  structure  in  the  upper  ocean.  Many  of  these 
data  have  been  made  possible  by  the  widespread  use  of  towed  vehicles  and  ship- 
board Doppler  sonars.  We  review  what  these  observations  have  revealed  about 
processes  in  the  upper  ocean,  with  special  emphasis  on  mid-latitudes  and  hori- 
zontal scales  of  1-100  km.  Fronts  are  one  of  the  most  prominent  features  of  the 
mixed  layer  on  these  scales.  The  vertical/across-ffont  circulation  created  at  a 
front  may  be  quantified  assuming  quasigeostrophic  dynamics.  Typical  of  the 
mixed  layer  are  fronts  that  are  warm  and  salty  on  one  side,  and  cool  and  fresh  on 
the  other  such  that  the  density  contrast  across  the  front  is  small.  The  ubiquitous 
existence  of  such  compensated  fronts  suggests  that  horizontal  mixing  in  the 
mixed  layer  is  an  increasing  function  of  the  horizontal  density  gradient.  Argu- 
ments based  on  the  conservation  of  potential  vorticity  indicate  that  1 0-km-scale 
(roughly  the  Rossby  radius  of  deformation)  eddies  in  the  mixed  layer  should  be 
predominantly  cyclonic.  Recent  observations  of  mixed-layer  vorticity  are  posi- 
tively skewed. 


Introduction 

The  database  of  small  horizontal  scale  upper  ocean 
observations  has  grown  significantly  in  recent  years  due  to 
the  increased  use  of  towed  vehicles  (such  as  SeaSoar), 
shipboard  Acoustic  Doppler  Current  Profiler  (ADCP)  and 
Global  Positioning  System  (GPS)  navigation.  A brief  re- 
view of  these  measurements  is  the  purpose  of  this  paper, 
with  special  emphasis  on  the  mid-latitude  open  ocean  on 
scales  of  1-100  km. 

Relevant  characteristics  of  the  instruments  used  in  upper 
ocean  surveying  are  as  follows.  SeaSoar  is  a towed  vehicle 
capable  of  cycling  from  the  surface  to  over  300  m every  3- 
4 km  at  a tow  speed  of  4 m s'1.  SeaSoar  usually  carries  a 
CTD,  providing  measurements  of  temperature,  salinity, 
and  pressure.  A 150-kHz  ADCP  in  typical  use  allows  hori- 
zontal resolution  of  velocity  quite  similar  to  that  of  a Sea- 
Soar, with  vertical  coverage  of  200-300  m.  High  quality 
GPS  navigation  is  a key  part  of  the  surveying  system,  es- 
pecially for  retrieving  the  best  absolute  velocity.  P-code 
receivers,  in  use  on  many  research  vessels,  offer  positional 
accuracy  of  about  3 m. 

Observations  are  used  to  address  three  issues.  (1)  The 
anatomy  of  an  upper  ocean  front  is  discussed,  including 
the  vertical/across-ffont  circulation  and  associated  heat 
flux.  (2)  Compensation  of  horizontal  temperature  and  sa- 
linity gradients  in  the  mixed  layer  is  shown  common.  (3) 


Cyclonic  vorticity  is  preferred  at  small  horizontal  scales  in 
the  mixed  layer.  For  each  of  the  above  issues,  we  discuss 
what  the  observations  suggest  about  mixing  in  the  ocean. 

The  anatomy  of  an  upper  ocean  front 

The  most  common  use  of  the  SeaSoar/ADCP/GPS  sys- 
tem has  been  to  survey  individual  features  in  the  upper 
ocean.  Such  a survey  of  the  Azores  Front  in  the  North  At- 
lantic, done  in  March  1992,  demonstrates  many  phenom- 
ena common  to  oceanic  fronts  [Rudnick  and  Luyten,  1996]. 
The  survey,  consisting  of  a series  of  north/south  sections 
separated  by  one-quarter  of  a degree  in  longitude,  reveals  a 
sharp  density  front  (Figure  1).  Horizontal  velocity  ap- 
proaches 0.5  m s'1,  and  is  generally  directed  parallel  to 
isopycnals.  This  intuitive  result  has  two  implications.  First, 
surveying  techniques  are  accurate  enough  that  coherent 
velocity  and  density  fields  can  be  measured.  Second,  and 
more  importantly,  the  rate  of  change  and  vertical  advection 
of  density  must  be  small.  The  velocity  field  is  in  the  form 
of  a jet  strongest  at  the  density  front.  The  jet,  and  its  asso- 
ciated straining,  causes  streamers  in  the  alongfront  direc- 
tion. Enhanced  horizontal  variability  is  typical  of  frontal 
regions  [ Ferrari  and  Rudnick,  2000]. 
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Figure  1.  Potential  density  and  horizontal  velocity  in  the  mixed 
layer  at  the  Azores  Front.  The  contour  interval  is  0.025  kg  m'\ 
and  the  velocity  scale  is  shown.  Note  the  jet  directed  along  the 
density  front. 
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The  vertical/across-front  structure  is  seen  in  a 
north/south  section  along  24.25°W  (Figure  2).  Mixed-layer 
depth  trends  from  about  100  in  deep  in  the  northern  part  of 
the  section  to  roughly  150  m in  the  south.  The  mixed-layer 
base  is  marked  as  a region  of  strong  stratification.  The 
front  is  apparent  as  isopycnals  slope  upward  and  to  the 
north  across  the  mixed  layer.  One  can  think  of  the  structure 
almost  as  a piece  of  the  mixed-layer  base  peeling  off  and 
intersecting  the  surface.  Because  the  survey  is  from  late 
winter,  following  a period  of  relative  calm,  restratification 
is  just  beginning.  A shallow  stratified  layer  lies  over  the 
mixed  layer  north  of  the  front.  Isopycnals  retain  a strong 
slope  in  the  thermocline  beneath  the  surface  front.  The 
Azores  Front  is,  in  fact,  a permanent  feature  that  penetrates 
to  a depth  of  2000  m [Gould,  1985].  There  is  some  sugges- 
tion of  subduction,  as  mixed-layer  water  from  the  north 
side  of  the  front  appears  to  slide  southward  and  below  the 
outcropping  front. 

The  simplest  dynamical  balance  one  can  imagine  testing 
is  thermal  wind.  Such  a test  is  possible  by  comparing  the 
independently  measured  fields  of  velocity  shear  and  den- 
sity (Figure  2).  Strong  shear  of  the  correct  sign  is  apparent 
along  the  outcropping  isopycnals  in  the  mixed  layer,  and 
downward  and  southward  beneath  the  mixed  layer.  Alter- 
nating bands  of  positive  and  negative  shear  are  apparent  at 
the  mixed  layer  base  north  of  the  front.  These  shears  are 


clearly  out  of  thermal  wind  balance,  and  are  likely  near- 
inertial. 


Latitude 


Figure  2.  North/south  section  of  potential  density  and  vertical 
shear  of  velocity.  The  contour  interval  of  density  is  0.025  kg  nv\ 
Note  the  high  shear  along  the  sloping  front. 


The  vertical  velocity  field  may  be  inferred  assuming 
quasi-geostrophic  dynamics.  This  assumption  is  reasonable 
because  the  front  is  observed  to  be  nearly  in  thermal  wind 
balance.  Derivation  of  the  resulting  omega  equation 
[Hoskins  el  at.,  1978]  relies  on  the  elimination  of  the  time 
derivative  terms  in  the  momentum  and  density  equations. 
The  omega  equation  can  be  written 

V2  (A^U’W/3  = 2V  Q , (I) 

0 2 


where 


g_ 

Po 


f <9u  du  x 


dx 


d v 
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and  u4,  is  geostrophic  velocity,  the  divergence  operator  is 
two-dimensional  and  other  symbols  have  their  usual  mean- 
ing. Forcing  of  the  omega  equation  is  a function  of  the 
density  and  geostrophic  velocity  fields. 

The  elliptic  omega  equation  is  solved  for  the  survey  of 
the  Azores  Front  [ Rudnick , 1996].  A band  of  downwelling 
exists  near  and  to  the  north  of  the  front  (Figure  3).  Upwell- 
ing  appears  on  the  south  side  of  the  front  especially  at 
about  24°W.  A region  of  strong  downwelling  occurs  in  the 
SW  comer  where  the  flow  is  southward.  The  tendency  for 
downwelling  on  the  dense  (cold)  side  of  the  front  and  up- 
welling  on  the  less  dense  (warm)  side  results  in  a net  heat 
flux  averaged  over  the  horizontal  area  of  the  survey.  The 
heat  flux  peaks  at  the  base  of  the  mixed  layer  at  about  15 
W m"2.  This  heat  flux  is  comparable  to,  or  larger  than,  cli- 
matological values  compiled  on  1°  grids  [Marshall  et  al., 
1993]. 
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Figure  3.  Vertical  velocity  and  potential  density  in  the  mixed 
layer  at  the  Azores  Front.  The  contour  interval  for  potential  den- 
sity is  0.025  kg  m'3.  The  white  contour  indicates  zero  vertical 
velocity.  Note  the  band  of  downwelling  in  and  to  the  north  of  the 
front. 

A simple  theory  of  a semi-geostrophic  front  [ Hoskins 
and  Bretherton,  1972]  helps  to  understand  the  observa- 
tions. The  two-dimensional  model  relies  on  the  assumption 
that  the  across-front  momentum  balance  is  exactly 
geostrophic,  but  the  alongfront  balance  includes  horizontal 
and  vertical  advection  (see  Appendix).  The  front  is  formed 
by  a large-scale  deformation  of  an  initially  vertically  uni- 
form layer  with  weak  horizontal  density  gradients.  The 
solution  (Figure  4)  is  shown  at  time  4/a,  where  a is  the 
rate  of  strain.  The  large-scale  strain  estimated  from  the 
velocity  field  of  Figure  1 is  approximately  10'6  s'1,  so  the 
dimensional  time  would  be  46  days.  Sloping  isopycnals 


Figure  4,  The  semi-geostrophic  model  of  frontogenesis  at  time 
4 la.  Black  lines  are  isopycnals  with  a contour  interval  of  0.2  kg 
m'3.  Lines  of  constant  along-front  velocity  (blue)  have  an  inter- 
val of  0.1  m s'1.  The  red  lines  show  the  positions  of  particles  ini- 
tially at  10  and  50  m. 
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are  pinched  at  the  surface  forming  a front.  An  essential 
feature  of  this  model  is  that  an  infinitely  sharp  front  is 
formed  in  a finite  time.  Alongfront  velocity  is  in  the  form 
of  a surface  intensified  jet  with  the  dense  (positive  vortic- 
ity)  side  sharper  than  the  light  (negative  vorticity)  side.  As 
the  solution  is  Lagrangian,  it  is  a simple  matter  to  deter- 
mine the  position  of  particles.  The  initial  10  and  50  m sur- 
faces have  been  stirred  to  be  multi-valued  as  a function  of 
horizontal  position.  Strong  localized  downwelling  occurs 
on  the  dense  side  of  the  front  with  weaker,  more  diffuse, 
upwelling  on  the  less  dense  side.  This  pattern  of  vertical 
velocity  is  consistent  with  the  observations. 

Temperature/salinity  compensation 

Mixed-layer  thermohaline  fronts  arc  often  observed  to 
have  a gradient  from  warm  and  salty  to  cool  and  fresh  such 
that  the  density  contrast  across  the  front  is  small.  This  phe- 
nomenon is  termed  compensation  because  temperature  and 
salinity  compensate  in  their  effects  on  density.  Compensa- 
tion has  been  known  for  some  time  at  large  and  rneso 
scales  [Roden,  1975;  Roden , 1984],  More  recently,  com- 
pensation has  been  shown  to  exist  at  horizontal  scales  as 
small  as  10  m [Ferrari  and  Rudnick , 2000;  Rudnick  and 
Ferrari,  1999]. 

The  existence  of  compensation  in  the  mixed  layer  is 
consistent  with  theory  proposed  by  Young  [1994]  and  Fer- 
rari and  Young  [1997].  Consider  a vertically  mixed  layer 
with  an  initially  random  distribution  of  temperature  and 
salinity.  Temperature  and  salinity  gradients  compensate  in 
their  effect  on  density  in  some  regions,  while  in  others 
density  gradients  exist.  Density  gradients  slump  due  to 
gravity,  creating  sloping  isopycnals.  Vertical  mixing  re- 
sults in  weakened  horizontal  density  gradients.  This  proc- 
ess is  essentially  shear  dispersion  [Taylor,  1953]  where  the 
shear  is  due  to  slumping  density  gradients,  and  the  mixing 
is  caused  by  any  of  the  processes  that  make  the  mixed 
layer  vertically  uniform.  Temperature/salinity  gradients 
that  are  compensated  do  not  slump  and  therefore  do  not 
experience  shear  dispersion.  The  net  effect  is  that  density 
gradients  diffuse  while  compensated  gradients  persist. 

Compensation  is  manifest  as  a strong  correlation  be- 
tween temperature  and  salinity  in  a horizontal  SeaSoar  tow 
in  the  N.  Pacific  winter  mixed  layer  [Ferrari  et  al.,  this 
volume].  Temperature  and  salinity  fluctuations  coincide  at 
all  resolved  scales.  The  observations  are  thus  consistent 
with  a horizontal  diffusivity  that  is  an  increasing  function 
of  horizontal  density  gradient. 

An  important  question  is  under  what  conditions  is  com- 
pensation typical  of  the  world’s  oceans.  This  question  is 
addressed  using  SeaSoar  data  from  the  midlatitude,  north- 
ern hemisphere  Pacific,  Atlantic,  and  Indian  Oceans  during 
all  seasons  [Rudnick  and  Martin,  2001].  The  relevant 
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quantity  is  the  density  ratio  R,  defined  as  the  ratio  of  the 
effect  of  a horizontal  change  of  temperature  on  density 
divided  by  that  of  salinity  on  density.  For  a compensated 
therniohaline  front,  R=  1. 


Figure  5.  The  conditional  probability  density  function  (pdf)  of 
mixed-layer  density  ratio  as  a function  of  mixed-layer  depth  from 
the  global  ocean.  Each  row  of  bars  represents  the  conditional  pdf 
for  the  given  mixed-layer  depth.  Height,  and  color  shading,  of  the 
bars  indicate  the  pdf  magnitude.  The  density  ratio  is  typically  I 
for  mixed  layers  deeper  than  75  m. 


The  conditional  probability  density  function  (pdf)  of 
mixed-layer  density  ratio  as  a function  of  mixed-layer 
depth  summarizes  results  from  the  global  ocean  (Figure  5). 
For  mixed-layers  deeper  than  75  m a density  ratio  near  1 is 
typical,  while  shallower  mixed  layers  have  a poorly  de- 
fined density  ratio.  The  shallowest  mixed  layers  are  tem- 
perature dominated  as  the  density  ratio  approaches  infinity. 
Mixed-layer  depth  may  be  considered  a proxy  for  the 
strength  of  vertical  mixing.  The  results  suggest  that  verti- 
cal mixing  is  essential  to  the  establishment  of  compensa- 
tion in  the  mixed  layer.  This  observation  is  consistent  with 
the  slumping  and  mixing  model  where  vertical  mixing  is 
an  essential  process.  On  the  other  hand,  compensation 
should  not  be  expected  in  shallow  mixed  layers  where  ver- 
tical mixing  is  weak. 

A density  ratio  of  two  is  typical  of  the  thermocline  be- 
neath the  mixed  layer  [Schmitt,  1981],  This  phenomenon  is 
believed  to  be  caused  by  salt  fingering  [Schmitt,  1994], 
What  is  remarkable  is  the  sharpness  of  the  transition  from 
R=1  to  R= 2 across  the  mixed-layer  base,  especially  when 
the  mixed  layer  is  deep  [Ferrari  and  Rudnick,  2000;  Rud- 
nick and  Martin,  2001].  The  change  in  density  ratio  from 
the  mixed  layer  to  the  thermocline  is  a reflection  of  the 
different  mixing  processes  working  in  the  two  regions. 


Mixed-layer  vorticity 


Rapidly  rotating  vortices  in  the  ocean  arc,  in  theory, 
predominantly  cyclonic.  A simple  explanation  invokes  the 
conservation  of  potential  vorticity  as  in  the  two- 
dimensional  Hoskins  and  Brcthcrton  scniigeostrophic 
model  (Figure  4).  The  potential  vorticity  for  this  model  is 


dp  d v 

YxTz  ' 


(3) 


The  initial  state  of  the  model  is  one  of  zero  q as  the  layer  is 
vertically  uniform  and  velocities  arc  zero.  As  the  front 
sharpens,  gradients  of  density  and  velocity  increase,  but  q 
must  remain  zero.  Because  of  thermal  wind,  the  second 
term  above  is  always  positive  (in  the  northern  hemisphere). 
To  conserve  q then,  the  first  term  must  be  negative.  As 
static  stability  requires  dp/dz<0,  absolute  vorticity  must 

remain  positive.  If  the  relative  vorticity  is  to  be  the  magni- 
tude of/or  larger,  it  must  thus  be  cyclonic. 

The  discussion  above  has  relevance  to  the  symmetric  in- 
stability [Hoskins,  1974],  The  symmetric  instability  con- 
siders an  exactly  circular  vortex,  so  a two-dimensional 
model  (radial/vcrtical)  is  sufficient.  Instability  occurs  if  the 
potential  vorticity,  given  by  (3),  is  positive.  Hoskins  points 
out  that  a statically  stable,  motionless  state  cannot  evolve 
frictional ly  and  adiabatically  to  a state  subject  to  symmet- 
ric instability. 

The  relevant  length  scale  in  the  ocean  is  the  Rossby  ra- 
dius of  deformation  for  the  mixed  layer,  of  order  10  km. 
The  impediment  to  making  accurate  measurements  of 
small-scale  vorticity  using  shipboard  ADCP  data  has  been 
navigation.  With  the  advent  of  military  grade  (P-codc) 
GPS,  this  problem  has  been  solved. 

Vorticity  is  estimated  using  data  of  3-krn  horizontal 
resolution  from  a cruise  in  the  North  Pacific  during  the 
winter  of  1997  [Rudnick,  2001].  Because  a SeaSoar  was 
used,  underway  estimates  of  mixed-layer  depth  are  avail- 
able. The  pdf  of  mixed-layer  relative  vorticity  normalized 
by  planetary  vorticity  has  a mode  near  zero  and  is  posi- 
tively skewed  (Figure  6).  The  skewness  is  apparent  as  a 
relatively  thick  tail  at  a normalized  vorticity  of  0.5.  The 
skewness  of  vorticity  is  0.36,  and  is  considered  statistically 
significant.  The  skewed  vorticity  is  likely  the  velocity  sig- 
nature of  spiral  eddies  seen  in  photographs  of  the  sea  sur- 
face [Munk  et  al.,  2000],  The  contribution  of  these  eddies 
to  mixing  in  the  ocean  is  an  open  question. 
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Figure  6.  Probability  density  function  (pdf)  of  relative  vorticity 
normalized  by  planetary  vorticity  in  the  mixed  layer.  A normal 
distribution  with  zero  mean  and  standard  deviation,  0. 1 8,  equal  to 
the  observations  is  shown  (solid  line).  The  skewness,  equal  to 
0.36,  is  apparent  as  the  relatively  high  value  of  the  pdf  near  a 
normalized  vorticity  of  0.5. 
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Here,  y is  the  along-front  direction,  z is  vertical,  p0  is  a 
constant  reference  density,  and /is  a constant  Coriolis  pa- 
rameter. Ertel’s  potential  vorticity  (3)  is  conserved  follow- 
ing fluid  particles 


Summary 


Assume  a solution  consisting  of  a geostrophic  confluent 
flow  plus  a two-dimensional  component, 


Shipboard  surveys  using  SeaSoar  and  ADCP  have 
greatly  improved  our  knowledge  of  small-scale  horizontal 
variability  in  the  upper  ocean.  Oceanic  fronts  are  observed 
to  be  in  approximate  thermal  wind  balance.  A restratifying 
across-front  circulation  can  be  inferred,  implying  a heat 
flux  of  order  1 0 W m'2  when  averaged  over  a 1 00  km2  re- 
gion surrounding  the  front.  Temperature/salinity  compen- 
sation is  typical  in  the  mixed  layer,  especially  when  the 
mixed  layer  is  deep.  The  observed  compensation  is  consis- 
tent with  the  slumping  and  mixing  model  of  Young  and 
colleagues,  and  a horizontal  diffusivity  that  is  a growing 
function  of  density  gradient.  Mixed-layer  vorticity  is 
skewed  positive  because  the  largest  vorticities  are  positive. 
The  observed  skewness  is  consistent  with  potential  vortic- 
ity conservation. 

Appendix 

The  following  is  a review  of  the  Hoskins  and  Bretherton 
[1972]  two-dimensional  semi-geostrophic  model  of  fron- 
togenesis.  A novelty  is  that  the  solution  is  in  fully  Lagran- 
gian  coordinates,  in  contrast  to  the  mixed  Eule- 
rian/Lagrangian  solution  in  Hoskins  and  Bretherton. 

The  semi-geostrophic  equations  appropriate  for  a front 
in  the  x-direction  are 


Li  - -ax  + u'(x,z,t ) (A7) 

v = ay+v(x,z,t ) (A8) 
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p = p(x,z,t) 
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(A10) 


w = \v{x,z,t)  ■ 


(All) 


Note  that  the  parameter  a may  be  a function  of  time. 
Substitution  of  (A7-1 1)  into  (Al-5)  yields  the  equations 
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du  dw  . . . 

1-  — = -a 

dx  dz 
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where 

M=v+fx. 

These  together  with  conservation  of  density  (A4)  are  the 
equations  we  need  to  solve. 


(Al) 
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The  prognostic  equations  (A  12)  and  (A4)  can  be  solved 
immediately  in  Lagrangian  coordinates,  with  no  explicit 
dependence  on  along-front  position, 

M = M0(x0,z0)e'p  (A  15) 

P = P(*o>2o)>  (A  16) 


where 


and  x0  and  z0  are  the  initial  across-front  and  vertical  posi- 
tions of  particles. 

A useful  identity  expresses  the  Eulerian  spatial  gradient 
of  a dependent  variable  h in  Lagrangian  coordinates. 


d(h,z) 

d(x,z) 

d(x0,z0) 

d(x0,z0) 

(A17) 


where  the  subscript  z on  the  lhs  makes  clear  that  z is  held 
constant.  Using  (A17),  (A13-14)  can  be  written  in  Lagran- 
gian coordinates  as 


d(x’z)  _.-g 
d (xo>zo) 


(A  18) 


x = z-S- 


dp 


f'Po  dx0 


ep  + A(x0,t), 


(A23) 


where  A(x0  ,t)  is  to  be  determined  from  boundary  condi- 
tions. Derivatives  of  x with  respect  to  x0  and  z0  can  be  de- 
rived from  (A23)  and  substituted  into  the  continuity  equa- 
tion (A1 8).  Integrating  the  result  with  respect  to  z0  gives 


dx0  2f2p0  dx] 


= z0e'p +B(x0,t). 


(A24) 


We  use  boundary  conditions  that  particles  on  the  top 
and  bottom  boundaries  remain  there 


z(z0  =±H)-±H . (A25) 


Evaluating  (A24)  at  the  boundaries,  and  using  (A23), 
yields  expressions  for  the  integration  constants 


A=x0e 


,-P 


B = H2--?  -d  -^ep . 
2f~Po  dxl 

The  solutions  for  x and  z are  thus 


x = xQe  +z 


f^Po  dx 0 


(A26) 

(A27) 


(A28) 


_£(x^  = _j__£(p£)_ 

d{x0,z0)  pod(x0,z0)' 


(A  19) 


z-z„  + 


2/:P0  dx( 


(A29) 


The  solution  in  Lagrangian  coordinates  (A  15-16)  is 
simple.  The  challenge  is  the  transformation  back  to  Eule- 
rian coordinates  (A18-19).  These  equations  are  nonlinear, 
and  possibly  difficult  to  solve  depending  on  the  initial  and 
boundary  conditions  used. 

Potential  vorticity  may  be  written 

(a20) 

Po  d(x,z)  p0  d(x0,z0) 

and  is  necessarily  conserved,  as  can  be  readily  seen  by  its 
expression  in  Lagrangian  coordinates. 

Consider  a layer  initially  vertically  uniform  in  density. 
Suppose  also  that  the  initial  horizontal  density  gradient  is 
so  small  that  along-front  velocity  and  its  vertical  shear  are 
negligible.  Using  these  initial  conditions,  solutions  are 

M = fxoe~0  (A21) 

p = p(x0).  (A22) 

Substituting  these  into  the  thermal  wind  equation  (A19) 
and  integrating  yields 


The  solution  in  Figure  4 has  an  initial  density  profile  that  is 
an  arctan  with  a length  scale  of  500  km,  in  a layer  200  m 
deep. 

A few  features  of  the  solution  are  worthy  of  note.  The 
absolute  vorticity  is 


+ z 


g d~P 

f2Po  dxl 


(A30) 


Vorticity  becomes  infinite  in  a finite  time,  provided  the 
initial  density  profile  has  nonzero  curvature.  This  will  hap- 
pen first  on  the  boundaries. 

The  slope  of  an  isopycnal. 


- g dP  .0 
dz  J,  f2p0dx  o 


(A31) 


depends  only  on  x0  and  t,  so  the  isopycnal  is  a straight  line. 
The  Richardson  number  is  given  by 


Ri  = 


N2 

(dv/dzy 


f 

dM/dx 


(A32) 


Therefore,  (A30)  implies  that  Ri  becomes  zero  in  a finite 
time  in  the  region  of  strong  vorticity  and  subduction  at  the 
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front.  Ri  is  1 elsewhere.  Mixing  must  become  important  as 
the  front  sharpens. 
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